NAVIGATIONAL SYSTEM AND METHOD UTILIZING 
SOURCES OF PULSED CELESTIAL RADIATION 

BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION 

The subject invention relates to a system and method for navigation 
utilizing sources of pulsed celestial radiation. In particular, the present invention 
directs itself to a mobile receiver for detecting pulsed signals generated by 
celestial sources of pulsed radiation. The mobile receiver is mounted on a 
spacecraft, satellite, planetary rover, or other vehicle and the received pulses are 
used for the calculation of navigational data for the vehicle. More particularly, 
this invention directs itself to a clock or timer in communication with the mobile 
receiver for generating a timing signal corresponding to the reception and 
detection of the pulsed celestial radiation. The timing signal is used to calculate 
a time offset between predicted and measured pulse reception at the mobile 
receiver. Incorporating predetermined models of the pulse arrival times within 
an inertial reference frame, a position offset of the mobile receiver can then be 
determined using the time offset. 

Further, the mobile receiver is in communication with a digital memory 
system. The digital memory system stores information including the positions 
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and pulse timing model parameters of known sources of pulsed celestial 
radiation with respect to the chosen inertial reference frame. Additionally, this 
invention directs itself to a navigational system including a processing means for 
calculating navigational data of the spacecraft, satellite, planetary rover, or 
vehicle based upon the calculated time offset in combination with the known 
positional and pulse timing model data of the sources of celestial radiation stored 
within the digital memory system. 
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PRIOR ART 

Celestial objects, or sources, have been utilized throughout history as 
navigational aids. The motions of the Sun, Moon, and planets, as observed from 
the Earth's surface, have provided the concept of time, and are a form of a 
celestial clock. Using catalogued almanac information, the observation of 
visible stars have provided travelers on the Earth a means to determine position 
information relative to observation stations fixed on the Earth. As these methods 
have matured over time, and with the addition of instrumented time clocks, or 
chronometers, the performance of navigation methods have improved. 

Navigation of vehicles above and beyond the Earth's surface have 
benefitted from this knowledge of celestial source-based navigation. Satellites 
orbiting the Earth and spacecraft traveling throughout the Solar System have 
relied on celestial sources to successfully complete their missions. Additionally, 
celestial source navigation systems have been augmented with human-made 
systems to further increase vehicle and spacecraft navigation performance. In 
fact, a wide variety of methods have been used to compute the navigation 
information of spacecraft which have traveled around the Earth, through the 
Solar System, and beyond the Solar System's outer planets, as far as the 
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heliopause. 

Navigation of spacecraft is defined as determining the spacecraft's three- 
dimensional position, velocity, and attitude at a specific time or times. Position 
determination of near-Earth missions have included the use of ground-based 
radar and optical tracking, Earth-limb horizon sensors, Sun sensors, and Global 
Positioning System (GPS) receivers. These sensors use knowledge and 
observations of celestial objects and phenomena to determine position relative to 
the Earth. The GPS system produces signals from multiple transmitting satellites 
that allow a receiver to determine its position from the ranges to the transmitting 
satellites. Position determination of spacecraft on interplanetary missions have 
utilized Doppler radar range measurements, standard radio telemetry, and target 
body imagery. Attitude, or orientation, of spacecraft has typically been 
determined using magnetometers, gyroscopes, star cameras, star trackers, Earth- 
limb horizon sensors and Sun sensors. Time has often been determined using a 
clock on-board the spacecraft, or through periodic computer resets from ground 
control stations. 

Radar range and tracking systems have been the predominant system for 
tracking and maintaining continuous orbit information of spacecraft. In order to 
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compute position of spacecraft, radar range systems compute the range and/or 
the angular orientation angles to the spacecraft relative to an observing station. 
This is achieved primarily through the reflection of signals transmitted from an 
Earth observing station by the space vehicle structure and measurement of the 
transmitted signal round-trip time. Although this system requires no active 
hardware on the spacecraft itself, it does require an extensive ground tracking 
system and careful analysis of the measured data against an electromagnetically 
noisy background environment. Using the best known position of the observing 
station and processing multiple radar measurements, the vehicle's orbit 
parame^rs can be computed. The position of the vehicle can be propagated 
analytically ahead in time using standard orbital mechanics along with known 
models of solar system object's gravitational potential field; as well as any 
known disturbance or perturbations effects, such as object body atmospheric 
drag. This propagated navigation solution is then compared to subsequent radar 
measurements and the navigation solution is corrected for any determined errors. 
This process continues until a satisfactory navigation solution converges which 
is within the expedition's required parameters. However, vehicle maneuvers or 
any unaccounted-for disturbances will affect the trajectory of the vehicle. 
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Without exact knowledge of these maneuver dynamics or disturbance effects, it 
is necessary for the propagation and radar measurement comparison to continue 
throughout the flight. 

As the spacecraft moves further away from Earth observation stations, the 
error in navigational data increases. To achieve the necessary range 
determination, the radar system requires knowledge of the observation station's 
position on the Earth to great accuracy. With sophisticated surveying systems, 
including GPS, such accuracy may be achieved. However, even with this precise 
knowledge, the position measurement can only be accurate to a finite angular 
accuracy. The transmitted radar beam, along with the reflected signal, travels in 
a cone of uncertainty. This uncertainty degrades the position knowledge in the 
transverse direction of the vehicle as a linear function of distance. As the vehicle 
gets more distant, any fixed angular uncertainty reduces the knowledge of 
vehicle position, especially in the two transverse axes relative to the range 
direction. Even utilizing interferometry, using the difference between multiple 
signals compared at two ranging stations, the angular uncertainty rapidly grows 
above acceptable limits. 

Alternatively, many deep space spacecraft have employed active 

6 



transmitters to be used for navigational purposes. The radial velocity is 
measured at a receiving station by measuring the Doppler shift in the frequency 
of the transmitted signal. The spacecraft essentially receives a "ping" from the 
Earth observation station and re-transmits the signal back to the Earth. Although 
improvements in the range measurement are made utilizing such system, 
transverse axis errors still exist, and this method has errors that also grow with 
distance. 

Optical tracking measurements for spacecraft position and orbit 
determination are completed in a similar fashion as radar tracking. Optical 
tracking uses the visible light reflected off a vehicle to determine its location. 
Some optical measurements require a photograph to be taken and the vehicle's 
position is calculated after analysis of the photograph and comparison to a fixed 
star background. Real-time measurements using such systems are typically not 
easily achieved. Additionally, optical measurements are limited by favorable 
weather and environmental conditions. Because most missions have been 
oriented around planetary observation, augmentation to the ranging navigation 
system can be made within the vicinity of the investigated planet. By taking 
video images of the planet and comparing to known planetary parameters (such 

7 



as diameter and position), the video images can determine position of the 
spacecraft relative to the planet. Since, often the objective is to orbit the planet, 
only relative positioning is needed for the final phase of the flight and not 
absolute position. 

To aid the position determination process, an accurate clock is a 
fundamental component to most spacecraft navigation systems. On-board clocks 
provide a reference for the vehicle to use as its own process timer and for 
comparison to other time systems. Atomic clocks provide high accuracy 
references and are typically accurate to within one part in 10 9 - 10 15 over a day. 
As calculated by Melbourne in "Navigation Between the Planets", Scientific 
American, Vol. 234, No. 6, June 1976, pp. 58-74, in order to track radio signals 
at accuracies of a few tenths of a meter, a clock with nanosecond accuracy over 
several hours is needed. This requires the clock to be stable within one part in 
10 13 . As early chronometers helped improve navigation over the Earth's ocean, 
more accurate chronometers will help navigation through the solar system. 

To gain increased autonomy in spacecraft operation, it is desirable to 
develop systems other than the operationally intensive human-controlled radar 
and optical position determination methods. Seeking methods that employ 
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celestial sources, which provide positioning capabilities and do not require labor- 
intensive operations, remain attractive. 

In order to navigate using celestial objects, precise knowledge of their 
positions relative to a defined reference frame at a selected time epoch is 
required. Catalogued ephemeris information of solar system objects provide this 
position information. The Sun, Moon, and planets all translate within the solar 
system in a reference frame viewed from the Earth. Since the orbits of the 
objects are Keplerian, and nearly circular, this translational motion is nearly 
periodic, repeating after a certain elapsed time span. It is exactly this periodicity 
that leads to the concept of time. The motion of these objects, or their angular 
displacement in their orbit, can be interpreted as a clock measuring time. 

Although seemingly "fixed" with respect to a frame on the Earth, the 
visible spectrum stars can also provide a measurement of time and, therefore, can 
be interpreted as a navigational reference. In this case, however, it is the Earth or 
spacecraft that provides the time measurement, by rotating or translating with 
respect to the fixed background of stars. The extremely large distances to the 
stars in the Milky Way galaxy and other galaxies essentially create the illusion 
that the stars are fixed. Just as the solar system rotates, however, so does the 
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Milky Way rotate and the Galaxy translates with respect to neighboring galaxies. 
Thus, objects are continually speeding away and towards the Earth at all times. 
However, the motion of the stars is so slow compared to many other 
measurements of time that this motion is perceived as fixed. 

In addition to providing a measure of time, just as humans "triangulate" 
their Earth position relative to identifiable landmarks, it is conceivable to use 
persistent star light as markers for triangulating position. Observing known stars 
allows a spacecraft to initially estimate its orientation and begin a process of 
determining its position relative to another object. However, due to the large 
number of visible stars, detecting specific stars can be time-consuming due to 
necessary almanac database searches. Also, since there is no method of 
determining "when" the visible light was sent from the stars and because the 
stars are located so far away, determining range information from an individual 
star to help "triangulate" a spacecraft's position is problematic. Only during the 
instance of occultation, or when a known celestial body passes in front of a 
selected star while it is being viewed, allows a dependable method of position 
determination directly from star light. 

However, individual stars that do have a uniquely identifiable signal, 
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whose signals are periodic, can be utilized directly as celestial sources for 
navigation purposes. 

Astronomical observations have revealed several classes of celestial 
objects that produce unique signals. A particularly unique and stable source is 
generated by pulsars. It is theorized that pulsars are rotating neutron stars. 
Neutron stars are formed when a class of stars collapse, and from conservation of 
angular momentum, as the stars become smaller, or more compact, they rotate 
faster. Neutron stars rotate with periods ranging from 1.5 ms to thousands of 
seconds. A unique aspect of this rotation is that for certain classes of pulsars, the 
rotation can be extremely stable. The most stable pulsars have stabilities on the 
order of 10" 14 when measured on time scales of a year, which is comparable to 
the best terrestrial atomic clocks. No two neutron stars have been formed in 
exactly the same manner, thus their periodic signatures are unique. Because 
pulsars provide signals that are unique, periodic and extremely stable, they can 
assist in navigation by providing a method to triangulate position from their 
signals. Pulsars can be observed in the radio, optical, X-ray, and gamma-ray 
ranges of the electromagnetic spectrum. 

Downs, in "Interplanetary Navigation Using Pulsating Radio Sources", 
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NASA Technical Reports, N74-34150, Oct. 1, 1974, pp. 1-12, presented a 
method of navigation for orbiting spacecraft based upon radio signals from a 
pulsar. However, both the radio and optical signatures from pulsars have 
limitations that reduce their effectiveness for spacecraft missions. In order to be 
effective, optical pulsar-based navigation systems would require a large aperture 
to collect sufficient photons, since few pulsars exhibit bright optical pulsations. 
The large number of visible sources requires precise pointing and significant 
processing to detect pulsars in the presence of bright neighboring objects. This 
is not attractive for vehicle design. At the radio frequencies which pulsars emit 
(~ 100MHz - few GHz) and with their faint emission, large antennae (likely 25m 
in diameter or larger) would be required, which would be impractical for most 
spacecraft. Also, neighboring celestial objects, including the Sun, Moon, close 
stars, and Jupiter, as well as distant objects, such as supernova remnants, radio 
galaxies, quasars, and the galactic diffuse emissions are broadband radio sources 
that could obscure the weak pulsar radio signals, as shown by Wallace in "Radio 
Stars, What They Are and The Prospects for Their Use in Navigational System", 
Journal of Navigation, Vol. 41, September 1988, pp. 358-374. The low radio 
signal intensity from pulsars would require long signal integration times for an 
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acceptable signal to noise ratio. 

Chester and Butman in "Navigation Using X-Ray Pulsars", NASA 
Technical Reports, N81-27129, June 15, 1981, pp. 22-25, propose the use of 
pulsars emitting in the X-ray band as a better option for navigation. Antennae on 
the order of 0.1m 2 could be used for X-ray detection, which is much more 
reasonable than a large radio antenna. Additionally, there are fewer X-ray 
sources to contend with and many have unique signatures, which do not get 
obscured by closer celestial objects. One complicating factor is that many X-ray 
sources are transient in nature. The transient sources are only detectable at 
irregular intervals as a result of a modulation in the accretion rate onto the 
celestial source. The "steady" sources (such as all of the rotation-powered 
pulsars) would typically be used for navigation. By cataloging pulsar positions 
and recording their signal periodicity and identifying parameters, as well as the 
data from other types of pulsed celestial radiation sources, a table of candidate 
stars can be created for use in navigation. These catalogs can then be stored in 
data memory format for use by algorithm processes onboard vehicles that detect 
pulsed radiation signals. Maintenance of these catalogs, and timely 
dissemination of data updates, is required for a high performance navigation 
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system. 
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SUMMARY OF THE INVENTION 
The present invention provides for a system and method for navigation 
utilizing sources of pulsed celestial radiation. The system includes a mobile 
receiver for detecting signals of pulsed celestial radiation. The mobile receiver is 
in communication with a clock or timer for generating a timing signal 
corresponding to the reception of the pulses. A processing means calculates a 
time offset between predicted and measured time of reception of the pulsed 
celestial signal with respect to a chosen inertial reference frame. The processing 
means is further in communication with a digital memory system that stores 
positional data of known sources of pulsed celestial radiation. Using the position 
data along with the calculated time offset, the processing means can calculate 
navigational data for the spacecraft, satellite, planetary rover, or other vehicle 
upon which the mobile receiver is mounted. 

It is a principal object of the subject system and method for navigation 
utilizing sources of pulsed celestial radiation to provide a navigational system 
utilizing a displaceable receiver means for detecting signals generated by a 
plurality of sources of pulsed celestial radiation. 

It is a further objective of the subject invention to provide a timer means 
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in communication with the displaceable receiver means for generating time of 
arrival signals corresponding to the time of detection of each signal. 

It is a further objective of the subject invention to provide a navigation 
system utilizing sources of pulsed celestial radiation including a processor means 
for computing navigational data based upon a calculated distance between the 
displaceable receiver means and a reference point in a chosen inertial reference 
frame. 

It is an important objective of the present invention to provide a 
navigation system including a digital memory system in communication with the 
processor means, the digital memory system having positional data and pulse 
profile characteristics of known sources of pulsed celestial radiation stored 
therein. 

It is a further important objective of the present invention to provide a 
system and method for navigation utilizing sources of pulsed celestial radiation 
where the processing means computes navigational data of a vehicle upon which 
the displaceable receiver means is mounted based upon the calculated time offset 
and the known positional data of the celestial signal sources. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a schematic view of the subject system for navigation utilizing 

sources of pulsed celestial radiation; 

Fig. 2 is a data plot illustrating the stability of several atomic clocks; 
Fig. 3 is a data plot illustrating the time stability of two millisecond 

pulsars; 

Fig. 4 is a data plot illustrating measured signal intensity generated by a 
sample pulsar; 

Fig. 5 is a data plot illustrating measured intensity generated by a second 
known pulsar; 

Fig. 6 is a schematic view of the navigational system mounted on a 
spacecraft; 

Fig. 7 is a schematic illustration of the geometrical calculations carried 
out by the navigational system; 

Fig. 8a is a side view of the navigational system; 

Fig. 8b is a top view of the navigational system; 

Fig. 8c is a bottom view of the navigational system; 

Fig. 9 is a flow chart for the method of navigation utilizing a source of 
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pulsed celestial radiation; 

Fig. 10 is a plot showing the measured pulse profile from a short 
observation of the Crab Pulsar; and, 

Fig. 1 1 is a schematic view of an alternative embodiment of the subject 
system for navigation utilizing sources of pulsed celestial radiation. 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 

I. Overview 

Referring now to Fig. 1, there is shown a navigational system 10 for 
calculating position data of a satellite 12 or other displaceable receiving platform 
based upon a signal 14 generated by a celestial body 16. Although shown in Fig. 
1 as a pulsar, the celestial body 1 6 can be any celestial body or source of a pulsed 
electromagnetic signal. 

II. Pulsed Celestial Sources 

Although Earth-based navigational systems, such as the Global 
Positioning System (GPS) provide navigational data including time, latitude, 
longitude, altitude, and attitude of a body on Earth, such systems have a limited 
scope of operations for vehicles moving and operating relatively far from the 
planet Earth. The system 10, however, utilizes celestial sources, such as pulsars, 
to provide a stable, predictable and unique signature signal in order to operate 
and navigate beyond the influence of the planet Earth. 
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Pulsars were initially discovered by Hewish, et al. in 1968 during a sky 
survey of scintillation phenomena due to inter-planetary plasma in the radio 
frequency range of approximately 100 MHz. Among the expected random 
noises emerged signals timed at regularly spaced intervals having periods of 
approximately one second. These periods were soon established to an accuracy 
of six or seven digits, making them one of the best determined astronomical 
constants outside the solar system. Although the exact nature of pulsars is not 
yet known, it is theorized that pulsars are rotating neutron stars. 

A neutron star is the remnant of a massive star that has exhausted its 
nuclear fuel and undergone a supernova explosion. When the remaining mass of 
the core is between 1 .4 and three times the mass of the sun, it collapses onto 
itself under its own gravity and a small extremely dense object is formed. This 
object lacks the necessary mass to become a black hole and it remains a stellar 
object formed largely of neutrons. Neutron stars are believed to be between 1 8 
and 30 km in diameter. Formed with an outer crust, the neutron stars are 
theorized to contain a neutron superfluid. Due to conservation of angular 
momentum, as these objects collapse, they begin to rotate at a faster rate. The 
young neutron stars typically rotate with periods on the order of tens of 
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milliseconds while older neutron stars eventually slow down having periods on 
the order of several seconds. A unique aspect of this rotation is that it can be 
extremely stable and predictable. The variance from the predicted spin rate can 
be as low as one part in 10 14 , which is comparable to the accuracy found in 
atomic clocks. 

Neutron stars further produce immense magnetic fields. Under the 
influence of these strong fields, charged particles are accelerated along the field 
lines to extremely high energies. As these charged particles move in the pulsar's 
strong magnetic field, powerful beams of electromagnetic waves are radiated out 
from the magnetic poles of the star. If the neutron star's spin axis is not aligned 
with its magnetic field axis, then an observer will sense a "pulse" of high energy 
photons as the magnetic pole sweeps across the observer's line of sight to the 
pulsar. It is theorized that these neutron stars and their respective pulses form 
what we know as pulsars. Since no two neutron stars are formed in exactly the 
same manner, the pulse frequency and shape produce a unique identifying 
signature for each pulsar. Due to their unique pulses, pulsars can act as natural 
beacons, or "celestial lighthouses", on an intergalactic scale, thus allowing for 
their use as navigational aids, such as in system 10 of Figure 1. 
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Many of the X-ray pulsars are "rotation-powered" pulsars, which are 
isolated neutron stars whose energy source is the stored rotational kinetic energy 
of the star itself. In addition to these rotation-powered pulsars, two other types 
of pulsars exist: "accretion-powered" pulsars and "anomalous pulsars". 
Accretion-powered pulsars are neutron stars in binary systems where material is 
being transferred from the companion star onto the neutron star. This flow of 
material is channeled by the magnetic field of the neutron star onto the poles of 
the star, which creates hot spots on the star's surface. The pulses are a result of 
the changing viewing angle of these hot spots as the neutron star rotates. These 
accreting pulsars are often subdivided into those with a high mass (typically 10- 
30 solar masses) or low mass (typically less than one solar mass) companions. 
The anomalous X-ray pulsars are powered by the decay of their immense 
magnetic fields (approximately 10 14 -10 15 Gauss). 

In addition to pulsars, such as pulsar 16 shown in Fig. 1, other variable 
pulsed celestial objects exist and can be used for spacecraft navigation. Many of 
these objects are highly variable in intensity and often employ different energy 
sources for their X-ray emissions. Such sources include active galactic nuclei, 
algol type stars, atoll sources, binary variable stars, black hole candidates, 
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bursters, coronal sources, cataclysmic variables, galaxy clusters, galactic ridge 
emissions, globular clusters, soft gamma repeaters, supernova remnants, white 
dwarf stars, and Z sources. Although these alternate X-ray sources have some 
desirable characteristics for use in navigation, some do not have detectable 
periodic or variable pulsations and others introduce additional complexities for 
signal processing and time-of-arrival determination. Whereas pulsars have 
reduced complexity, and are thus preferred as navigation beacon sources, all 
sources that produce predictable variable patterns can be implemented into the 
spacecraft navigation scheme. 

Although nearly all rotation-powered pulsars are constant in intensity, the 
accreting pulsars and most other X-ray source classes often exhibit highly 
aperiodic variability in intensity that limits their usefulness for precise time and 
position determination. Additionally, many are unsteady, or "transient", sources 
with a large range of duty cycles. Sources with X-ray flares of significant 
magnitude are also occasionally detected, which are high intensity signals lasting 
for short periods. Since neutron stars contain a solid crust and a super fluid 
interior, exchanges of angular momentum between the two materials can cause 
unpredictable "star-quakes", or "glitches", which can significantly vary the spin 
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rates of the stars. A navigational system that utilizes pulsars would have to 
address the transient, flaring, and glitch aspects of these sources. Additionally, 
as the diffuse X-ray background would be present in all observations, a 
navigation system that utilizes pulsed X-ray celestial sources must address the 
presence of this background noise in its data processing. 

Although many pulsars radiate at radio, optical, X-ray and gamma-ray 
wavelengths, not all have this feature. As X-ray and gamma-rays are difficult to 
detect on the ground due to the absorption of these wavelengths by Earth's 
atmosphere, detection within these bands requires observation to be made above 
the atmosphere. This complicates the detection and parameter verification of X- 
ray sources. Additionally, many X-ray sources are faint and require sensitive 
instruments to detect. 

Unlike human-made systems, such as the GPS, the distances to X-ray 
sources are not known to an accuracy that would allow them to be utilized as a 
ranging mechanism. However, the Galactic latitude and longitude of the X-ray 
sources can be determined with high precision. Most sources are located within 
our Galaxy, but still are very far from our solar system. Also, with many sources 
located within the Milky Way galactic plane, there are a limited number of bright 
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sources that could provide off-axis triangulation for position determination. 

III. X-Ray Source Observations 

Since the discovery of the first non-solar cosmic X-ray source, numerous 
balloon, rocket, and satellite borne instruments have surveyed the sky. Table 1 
provides a list of X-ray source survey missions. Different X-ray missions have 
observed the X-ray sky in various energy ranges, depending on instrument 
characteristics or mission goals. 
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Mission Name 


Mission Operation 


Energy Range (keV) 


Vela 5B 


1969-1979 


3-750 


UHURU 


1970-1973 


2-20 


OSO-7 


1971-1974 


1-10,000 


Skylab 


1973-1979 


0.1-0.3 


ANS 


1974-1977 


0.1-30 


Ariel V 


1974-1980 


03.-40 


OSO-8 


1975-1978 


0.15-1000 I; 


SAS-3 


1975-1979 


0.1-60 


HEAO 


1977-1979 


0.2-10,000 j 


Einstein 


1978-1981 


0.2-20 


Hakucho 


1979-1985 


0.1-100 


Tenma 


1983-1984 


0.1-60 


EXOSAT 


1983-1986 


0.05-20 


Ginga 


1987-1991 


1-400 


Granat 


1989-1998 


2-100,000 


ROSAT 


1990-1999 


0.1-2.5 


ASCA 


1993-2001 


0.4-10 i 


DXS 


1993 


0.15-0.28 


RXTE 


1995-Present j 


2-250 


BeppoSAX 


1996-2002 


0.1-300 


USA 


1999-2000 


1-15 


XMM 


1999-Present 


0.1-15 


[1 Chandra 


1999-Present 


0.1-10 



Table 1 
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During the late 1970's, the HEAO mission detected 842 X-ray sources 
within the 0.2-10 keV range. The German X-ray observatory ROSAT completed 
the latest comprehensive all-sky survey of the X-ray sky in the year 2000. This 
mission has detected 18,806 bright sources (above 0.05 X-ray photon counts per 
second in the 0.1-2.4 keV range), and detected a significant number of sources 
(105,924 objects) in its faint all-sky X-ray survey. 

The Australia Telescope National Facility has recently completed the most 
comprehensive radio pulsar study to date in their Parkes Multibeam Pulsar 
Survey. This survey has increased the number of known radio pulsars from 558 
to over 1,400. Many of these newly detected radio pulsars most likely radiate in 
the X-ray band, which adds to the number of potential pulsar navigation sources. 
By comparing various catalogs and published works, celestial sources can be 
investigated for their potential use in navigation systems. 

X-ray sources exist throughout the sky. Distances to X-ray objects in our 
Galaxy range from several parsecs to thousands of parsecs (1 parsec = 3.262 
light years = 3.086 x 10 16 m.). Most of the rotation-powered pulsars are detected 
within the Galaxy, with the rest being in the Magellenic Clouds. Magnetic fields 
of these objects range from several thousand to 10 14 Gauss (by comparison, the 
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Sun's magnetic field is a few tens Gauss in magnitude). Pulse periods range 
from between 0.0016 to 0.53 seconds for the rotation-powered pulsars, and from 
0.034 to 10,000 seconds for accretion-powered pulsars. 

As some pulsars have been observed for several decades, it has been 
shown that the stability of their spin rates compares well to the quality of today's 
atomic clocks. Figs. 2 and 3 provide comparison plots of the stability of atomic 
clocks and several pulsars measured in a third-order Allan Variance. The data 
for the plot of Fig. 2 comes from D.N. Matsakis, J.H. Taylor, and T. Marshall 
Eubanks, "A Statistic for Describing Pulsar and Clock Stabilities", Astronomy 
and Astrophysics, Vol. 326, 1997, pg. 924-928. The data for the plot in Fig. 3 
comes from V.M. Kaspi, J.H. Taylor, and M.F. Ryba, "High-Precision Timing of 
Millisecond Pulsars. III. Long-Term Monitoring of PSRs B 1855+09 and 
&\931+2Y\AstrophysicalJournah Vol. 428, June 1994, pg. 713-728. 

An accurate clock is a fundamental component in a spacecraft 
navigational system. On-board clocks provide a reference for a vehicle's process 
timer and for comparison to other time systems, and are critical for on-board 
systems such as communications. Atomic clocks provide high accuracy 
references and are typically good to within one part in 10 9 -10 15 over a day. In 
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order to track the motion of radio signals at accuracies of a few tenths of a meter, 
a clock with nanosecond accuracy over several hours is needed. This requires 
the clock to be stable to within one part in 10 13 . Fig. 3 illustrates that several 
pulsars meet the stability requirement. Older pulsars, particularly those that have 
undergone a long period of accretion in a binary system that spins them up to a 
millisecond period (the so-called "recycled pulsars") have extremely stable and 
predictable rotation rates. 

IV. Pulse Profiles and Arrival Times 

Source pulse profiles vary in amplitude, duration, number of peaks, and 
stability depending upon the nature of the pulsar. This uniqueness can aid in the 
identification of an individual source. Typically, standard profile templates are 
created by observing a source over long durations and "folding", or 
synchronously averaging, at the pulse period. This folding process produces a 
pulse profile with a very high signal-to-noise ratio, and using this standard 
profile, characteristics of the pulse can then be determined, such as period length, 
amplitude and variability. 

Fig. 4 illustrates an image of a standard template for the Crab Pulsar (PSR 
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B0531+21) in the X-ray band (2-10 keV). The intensity of the profile is a ratio 
of count rate relative to average count rate. This image shows two cycles of the 
pulse, with one main pulse and a secondary pulse at a lower amplitude. The 
period of this pulse is approximately 33.4 milliseconds (Epoch 48743.0 MJD) 
and several period derivatives have been detected. Fig. 5 shows an image of the 
pulse profile of PSRB 1509-58. Two cycles are also shown for this pulsar, with a 
pulse period of approximately 150.23 milliseconds (Epoch 48355.0 MJD). 

Pulsars are timed by comparing measured pulse time of arrival (TO A) 
with those predicted by a model. A TO A measurement is made by observing a 
pulsar for a short period of time and producing a folded profile. The observed 
profile will differ from the standard profile template, s(t), by several factors. 
Typically, the observed pulse will vary by a bias, b 9 a scale factor, k, and some 
random noise T}(t), which in the case of X-ray observations is typically dominated 
by Poisson counting statistics. The relationship between the observed pulse p(t) 
and the standard template is given by: 

p(t) =b +k-s(t-r)+r](t) (1) 
where r is the time offset between the two profiles that are to be measured. The 
TO A determination process consists of Fourier transforming the observed and 
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standard profiles and minimizing the fitting values to compute the bias, scale 
factor, and TO A. The TO A determination process is described in further detail 
below. Additionally, depending on the chosen reference time for a specific 
navigation system application, the time measured by the clock 24 may need to be 
referenced to a terrestrial time standard, such as Coordinated Universal Time 
(UTC), or Temps Atomique International (TAI International Atomic Time) and 
the necessary time corrections are added accordingly. 

V. Measurement of TO A 

The fundamental measurable quantity of a pulsar-based navigation system 
is the time of arrival (TO A) of a pulse at the detector 22 as referenced to the 
clock 24. The detector 22 records the time of each individual X-ray photon with 
respect to the clock 24 to some high precision, on the order of one microsecond 
or better. Let the recorded time of the nth photon be t„. During a total 
integration time, T, a large number of photons N will have their arrival time 
measured by the clock 24 and recorded. Individual photon arrival times from % 
to t N are then converted to their equivalent time as the SSBC, as described in 
Section VII below, and "folded" at the predicted pulse period of the pulsar based 
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on the timing model for the pulsar. A binned pulse profile is then constructed by 
dividing the pulse phase into M equal bins and dropping each of the N photons 
into the appropriate phase bin. Such a measured profile is illustrated in Fig. 10. 

The TOA is then determined by measuring the offset between the phase of 
the measured profile with respect to the high signal-to-noise standard profile 
template for that source. This is based on the assumption that, after averaging a 
sufficiently large number of pulses, a pulse profile recorded in the same energy 
range is invariant with time. Thus, the phase zero point on the template profile 
sets the definition of phase zero for that source. The template can be aligned 
with an arbitrary point in the profile as phase zero, but two conventions are 
commonly used. Either the peak of the main pulse can be aligned as phase zero, 
or the profile can be aligned such that the phase of the fundamental component 
of its Fourier transform is zero. The latter method is preferred (which reduces to 
the former in the case of a single-symmetric pulse profile) because it is more 
precise and generally applicable, and allows for easy construction of standard 
templates by measuring the phase of the fundamental Fourier component and 
then applying a fractional phase shift to the profile. 

It is important to determine the TOA with an accuracy that is determined 
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by the signal-to-noise ratio of the profile, and not by the choice of the bin size. 
A standard cross-correlation analysis does not allow this to be easily achieved. 
However, the method given by J.H. Taylor, "Pulsar Timing and Relativistic 
Gravity", Philosophical Transaction of Royal Society (London A), Vol. 341, 
1992, pp. 1 17-134, is independent of bin size and can be implemented into the 
navigation system. The technique employs the time shifting property of Fourier 
transform pairs. 

The Fourier transform of a function shifted by an amount x is the Fourier 
transform of the original function multiplied by a phase factor of e 2**ir . Since 

the measured profile differs from the template by a time shift and a scale factor 
plus random noise, as in Eq. (1), it is straightforward to transform both the 
profile and the template into the Fourier domain. The parameters in Eq. (1) are 
then easily determined by a standard least-squares fitting method. The final 
measured TO A of the pulse is then determined by adding the fitted offset x to the 
recorded start time of the integration, to. 
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VI. Pulse Phase Timing Models 

In order to predict when pulses are expected to arrive at an observing 
station, a pulsar phase timing model must be developed. In order to compute 
accurate arrival times of pulses, measurements must be made relative to an 
inertial frame; i.e., a frame unaccelerated with respect to the pulsars. As most 
observations are typically made on the Earth, or in a spacecraft moving about the 
Earth, such as satellite 12 shown in Fig. 1, this data collected while in a moving 
frame must be transferred to an inertial frame. General relativity theory provides 
a means of transferring this data to an inertial frame, and these equations relate 
the emission time of photons that emanate from a source as the photons travel 
through curved space-time to the arrival at a station. 

The parameterized post-Newtonian (PPN) frame provides one option of 
an inertial frame and a reference time for accurate pulse timing. Several 
astronomical observatories use the PPN frame. Additionally, a frame exists that 
simplifies some of the general relativistic equations. This frame, the Temps 
Dynamique Baricentrique (TDB-Barycentric Dynamical Time), has its origin at 
the solar system barycenter (SSBC). This origin is just outside the surface of the 
Sun, as shown in Fig. 1, and is the center of mass of the solar system. Fig. 1 
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shows the relationships of pulses 14 from pulsars to the inertial TDB frame and 
to a spacecraft 12 orbiting the Earth 20. The positions of the spacecraft 12 and 
the center of the Earth 20 relative to the SSBC are shown. Earth-based 
telescopes or Earth-orbiting spacecraft can initially use terrestrial time standards 
to reference their observation. Then a transformation, derived below, can be 
utilized to convert from terrestrial time to TDB time. 

The phase, ^ , of arriving pulses can be detected at the SSBC using a 

pulsar-phase model given by 



known as the "pulsar spin equation", or "pulsar spin down law", where t is the 
time of observation, T 0 is the reference epoch, and co, a?, cb are the pulse 
frequency and its derivatives. This equation can also be represented using 
period, P, (also angular velocity Q = 2 mo) using the following relationships, 



*0 = <fi(T 0 ) + a(t- T 0 )+^cb(t- T 0 ) 2 + ^a>(t-T o y (2) 
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to produce 



*0 = XTj+jtt-Tj-^rit- T o y + ^_^L)(,_ t;) 3 (4) 

Eq. (4) allows the determination of the phase of a pulse signal at a future time t, 
relative to a reference epoch T 0 . Thus, at time it is possible to predict when the 
peak amplitude of the pulsar signal is expected to arrive at the SSBC. However, 
as no detectors exist at the SSBC, observation corrections must be made to allow 
Eq. (4) to be used for predicting pulsar pulse arrivals at the actual detector 
location. The model shown in Eq. (4) utilizes periods and its first two 
derivatives, however, any number of derivatives can be used depending on what 
is required to model a particular pulsar's timing behavior. Many more 
parameters may be required for some pulsars, such as those in binary systems. 
The accuracy of the model prediction depends on the quality of the known 
timing model parameters and on the intrinsic noise of the pulsar rotation. 
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VII. Reference Origin Offset Correction 

In order to adequately model pulsar signals, an inertial frame must be 
chosen for signal timing, and the SSBC can be utilized. A navigation system on 
a spacecraft, such as spacecraft 12, would be comprised of a sensor that would 
detect pulsar signals at the spacecraft location and a clock that would time these 
signals relative to a terrestrial time standard. Fig. 6 illustrates spacecraft 12 of 
Fig. 1 receiving the pulses 14 generated by pulsar 16. A pulse sensor 22 is 
mounted on spacecraft 12, as is a clock 24. If the spacecraft position and 
velocity are known, then pulse arrival time at the spacecraft can be transferred to 
the SSBC and compared to Eq. (4). 

From Fig. 1, it can be seen that from a spacecraft's position r relative to 
the SSBC at time f, the offset of time a pulsar signal arrives at a spacecraft 
compared to its arrival at the SSBC to first order is, 

- 

t b -t obs = At = — (5) 
c 

where t b is the time of pulse arrival at the SSBC, the t obs is the time of 
observation at the spacecraft, h is the unit direction from the origin to the pulsar, 
and c is the speed of light. Since many pulsars are so distant from the Earth, the 
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unit direction to the pulsars can be considered constant throughout the solar 
system. Closer pulsars must include the pulsar's motion effects and user position 
when determining their unit direction. The time used in Eq. (5) is referred to as 
coordinate time, or the time measured by a standard clock at rest in the inertial 
frame. A spacecraft's clock, unless actually at rest (zero- velocity) with respect to 
the SSBC and at the same gravitational potential, does not measure coordinate 
time. A spacecraft's clock measures proper time, or the time a clock measures as 
it travels along a four-dimensional space-time path. The goal of a pulsar-based 
navigation system would be in part to provide accurate position information of 
the spacecraft. This could only be accomplished by accurately timing pulsar 
signals and then correctly transferring this time to the SSBC. If a performance 
goal of such a navigation system is to provide accurate position information on 
the order the less than 300 meters, then the system must accurately time pulses to 
at least less than 1 \is (=300/c). To achieve these types of accuracies, general 
and special relativistic effects on a clock in motion relative to an inertial frame 
and within a gravitational potential field must be considered. A derivation of 
these effects on time is provided below. 

The four dimensions of the space-time coordinate frame can be 
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generalized to 



{ct,x,y,z} = {x°,x\x 2 9 x 3 } 



(6) 




with ct representing the dimension related to time, and |x,>%z} representing the 



spatial coordinates. A scalar metric, ds 2 , can be created which is invariant with 
respect to arbitrary transformations of coordinate. This can be defined as 



where g afi = g afi (ct 9 x 9 y,z) is a function of time and spatial coordinates and 

the elements of g a/3 form a symmetric, covariant tensor, and the dx* terms define 

the space-time path of an object. The path taken by a light ray or particle in 
space-time is called a world line. The proper time measured by a clock, r , as 
the clock moves along a world line in space-time frame, is related to the scalar 
invariant via 



ds 2 = g a0 dx a dx p 



(7) 



\ds\ 



(8a) 



or 
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(8b) 



From special relativity, the metric tensor can be written in Cartesian, flat- 
space (absence of gravity) as 

ds 2 = - c 2 dt 2 + (dx 2 + dy 2 + dz 2 ) (9) 

From general relativity in a weak-gravitational field and nearly flat space, this 
metric has been shown of the order 0(l/c) to be, 



ds 2 = -[\-^jc 2 dt 2 + (dx 2 + dy 2 + dz 2 ) 



(10) 



where U is the total gravitational potential acting on the spacecraft clock. Using 
the spatial velocity term of 



v 2 = 
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(11) 



which represents the speed of the spacecraft's local frame through the solar 
system, Eq. (10) can be divided by dt 2 to yield, 



ds 
dt 
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Using Eq. (8b) and taking the square root, Eq. (12) becomes, 
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(13) 



Using Taylor-series expansion on Eq. (13) and keeping only the terms of order 
0(l/c 2 ) yields, 



dr = 



1 c 2 2\cj 



dt 



(14) 



with a maximum error of 10" 12 s. 

By integrating Eq. (14) a solution of the proper time relative to the 
coordinate time can be determined for a spacecraft clock. 

Integrating Eq. (14) over time yields, 



dt (15) 



For an earth-orbiting spacecraft, the velocity of the spacecraft's clock can be 
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related to the earth coordinate frame to express v as, 



(16) 



where v E is the velocity of the Earth in the SSBC frame, an p sc/E is the 

spacecraft's velocity with respect to the Earth in the SSBC frame. Expanding 
Eq. (16) and ignoring small terms, Eq. (15) can be rewritten as, 



where v E =\\v E \\ . The third term on the right-hand side is often referred to as the 
Sagnac effect. Using the following representation of pulsar position: 



where D 0 is assumed to be the position at a fudicial time, T 0 , we can arrive at a 

solution for the second term on the right-hand side of Eq. (17). This model of 
pulsar motion assumes a constant velocity, V 0 , referred to as "proper-motion" of 

the pulsar, which includes radial and transverse speed terms. Using the 
additional assumptions that pulsars are very far away from the spacecraft such 



(*-/.)= (r-f„) + £ 




D=D 0 +V 0 (T-T 0 ) 



(18) 
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that D 0 » V 0 {T- T 0 ) , and that the direction to the pulsar is n « D Q I \\D 0 \\ , the 
solution becomes, 



{t - lc ) = (r-r a ) + (T-T e )-^(n-V 0 )A,-(n.r)}- ' [|f| J 
+ 0 + ^)Z ^[ ln |(" ' 0 ) + |fc I] + ;j( V Ac/* ) 



where r } is the position of the spacecraft relative to the jth planetary body in the 
solar system, and y is the time dilation factor. Since the values of V 0 are small, 

the Sun imposes the primary gravitional field within the solar system, and y can 
be set equal to one, the expression Eq. (19) can be further simplified as, 
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(t-t 0 )-(T-T 0 ) = (r-T 0 ) + 



nr 



c 2d A 



Win 12 



-(«r) 2 ] 



(20) 



where fi=GM s is the gravitational parameter of the Sun. Setting the left-hand 
side to be the change in time measured at the SSBC, Eq. (20) may be written as, 

= Ar + ^-^prj[H 2 - ("^ 2 )]+ ^H» r) + pr(v £ • Ac,*) ( 21 ) 

The second term on the right-hand side of Eq. (21) is the first order Doppler 
delay, and the third term is due to the effects of annual parallax. Together these 
two terms are referred to as Roemer delay. The fourth term is the Sun's Shapiro 
delay, and the fifth term is the relativistic effects of a portable clock in motion 
relative to the SSBC. With the addition of the fifth term, this equation matches 
well with those used by Earth-based pulsar timing observations. The interstellar 
medium dispersion measure term, appearing as a correction for radio 
observations, is considered zero in Eq. (21) for high frequency X-ray radiation. 
Given a measurement of a pulse arrival by the spacecraft's clock, the 
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measured time, r f must have its own reference time (for example, UTC or TAI) 
converted to TDB, then Eq. (21) can be used to determine the arrival time of that 
same pulse at the SSBC. Eq. (4) can be further utilized for a comparison of 
predicted and detected arrival times. Eq. (21) requires accurate ephemeris 
information to provide the SSBC location, the Sun's gravitational parameters, 
and the velocity of the Earth. It is intended for spacecraft in orbit about the 
Earth, and extensions to interplanetary spacecraft can be made through a similar 
derivation. 

VIII. Time Determination and Correction Method 

If accurate spacecraft position is known, Eq. (21) can be used to 
determine the expected arrival time of the pulse at the SSBC. If a spacecraft 
clock is in error, then the offset of the measured arrival time to the expected 
arrival time provides a measure of this clock error. 

Assume a spacecraft clock can be represented by the following, 

r T = r c + b+k(r c - r 0 )+ ^j(r c - T 0 f + t{t c ) (22) 
where r T is the true time, r c is the measured time from the clock, r 0 is a 
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reference time, b the clock bias, k is the clock rate, j is the clock jitter, and rj is 
the noise within the clock. Given estimates of b, k, and j 9 a Kalman filter can be 
created using this clock model. The measurement to this filter would be the 
offset between the estimated clock error and the computed clock error, 



where r P is the pulse arrival time from Eqs. (4) and (21). The frequency of the 
pulsar signal is used to stabilize a local clock using this filter. 

IX. Celestial Object Occulation Position Determination Method 

In the Earth-Limb occultation method, a detector is pointed towards the 
Earth's limb. As the vehicle revolves around the Earth in it's orbit, X-ray 
sources move behind the limb and then reappear on the other side. The time 
spent behind the Earth represents a chord length of the Earth's disk. Knowing 
the source position and the Earth's dimensions, it is possible to determine the 
position of the vehicle relative to the Earth. Knowledge of the Earth's 
atmosphere is required, since the X-ray signals would begin to be absorbed by 




(23) 



(24) 
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the atmosphere as the source passes close to the limb. 

Although this occultation method is described using the Earth's limb, any 
celestial body can be utilized in this concept. Bodies with known dimensional 
characteristics and positional ephemeredes are good candidates for this method, 
as long as they occult the path of a viewed X-ray source. Although we refer to 
X-ray sources, it is well understood that sources of any celestial radiation may be 
used. 



X. Position Estimator Method 

Due to the unique signatures of pulsars, it is possible to determine the 
position of a spacecraft, such as spacecraft 12 in Fig. 1. The position is 
determined relative to a desired inertial reference frame. Although the SSBC, 
shown in Fig. 1, provides one such frame, it is often more useful for mission 
operations to also relate vehicle position to the Earth's position. The occultation 
method described above and the position estimator and cold-start processes 
described below are methods of position determination relative to the Earth. 
Methods for determining the position of a spacecraft on interplanetary missions 
can be extended from these Earth-based examples. 
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In a "position estimator" method, pulsar signals received at a spacecraft 
are offset from those arriving at the SSBC primarily by the distance between the 
SSBC and the spacecraft, as in Eq. (5). Conversely, if the time offset can be 
directly measured, then the position offset can be computed. A pulsar's 
individual pulses arrive at the SSBC as per the model defined in Eq. (4). This 
equation allows the determination of when the Ath pulse will arrive at the SSBC 
relative to a chosen reference T 0 . As it moves away from the SSBC, a spacecraft 
sensor will detect a pulse at a time relative to T 0 based upon Eq. (21). However, 
if the spacecraft position is not known, then Eq. (21) cannot be used to transform 
spacecraft time to the SSBC. Pulses can still be detected and timed, but only 
relative to the spacecraft itself. 

In the position estimator scheme, a previously estimated value of 
spacecraft position is utilized. From this estimated position, r , the predicted 
arrival times at the spacecraft are compared to the detected arrival times. As 
shown in Fig. 7, the discrepancy in these values provides an estimate of the 
offset position, Sr . As shown in Fig. 7, the error in position will relate to a 
measured time offset of a pulse, in the line-of-sight to the pulsar. Using pulsars 
at different locations provides line-of-sight measurements in each pulsar's 

48 



direction. Combining these different pulsar measurements solves the position 
offset in three dimensions. 

Eq. (21) can be used to determine a position offset. Defining the offset 
Sr such that, 



r ~ r - Sr 



a time error can be computed by, 

*= t b (r 9 r 9 n 9 D 09 v E )- 7 b (r 9 r 9 n 9 D Q9 v E ) 



(25) 



(26) 



Approximations to the non- linear terms in Eq. (21) can be generated using Eq. 
(25) as, 



llFll 2 *\\r\\ 2 +2r-& 



(27) 
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(28) 
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(29) 



With the above approximations in Eqs. (27), (28), and (29), the expression in Eq. 
(21) can be placed in linear form with respect to Sr as, 
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(30) 



Eq. (30) can be rewritten in vector form as, 



(31) 



where the vector N = n(t , £) 0 , v £ ) is composed of the terms from Eq. (30). 



This linear expression can be assembled for k different pulsars, with sufficiently 
different locations to create a matrix of observations, 



(32) 
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The left-hand side of Eq. (32) is computed using predicted and actual arrival 
times of a pulse at the SSBC and on the spacecraft. Given a current estimate of 
r at?, Eq. (21) can be used to determine the current time at the SSBC, 7 b . 

Using Eq. (4), the expected arrival time of the maximum amplitude of phase of 
the next pulse can be determined, t b . The difference of these can then be 
computed, St = t b -7 b = expected arrival time - current time. On the spacecraft, 

the proper time of a pulse arrival can be measured, r , and the difference 
between actual arrival and current time can be calculated Sr = t-t . Using 
these computations, and the vectors N , Eq. (32) can be solved for Sr , and the 
estimated position of the spacecraft can be updated using Eq. (25). If no errors 
in modeling or measurement are present, then Eq. (32) solves for accurate 
position offsets. However, errors do remain which limit the performance of this 
system; 6t contains errors in pulsar modeling, Sr contains system level timing 
errors and pulse signal timing errors, and these errors along with any errors in 
parameters of N contribute to the errors in Sr , including pulsar position 
uncertainty and earth ephemeris accuracy. A phase cycle ambiguity is still 
present, as these equations only can relate to a fraction of a cycle and not 
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"which" cycle is being detected. Additionally, vehicle motion that is significant 
during the time span between current time and pulse arrival time must be 
addressed in an implementation of this position estimate system. A Kalman filter 
incorporating vehicle dynamics and a measurement model from Eq. (32) can be 
used to successfully update vehicle position estimates. Additional complexity is 
added if one chooses to incorporate binary pulsar observations, and these extra 
terms must be added to Eq. (32). 

XI. Cold-Start Position Determination Method 

Although determining an offset to an estimated position is useful for some 
systems, navigation systems that can determine absolute position using no prior 
knowledge are very desirable. To determine absolute position from the signals 
of pulsar 16, shown in Fig. 1, it is necessary to determine which integer signal 
cycle, or pulse period, is being detected at a certain time. Once it is determined 
which specific cycle is detected, then range between the SSBC and the spacecraft 
can be determined along that pulsar's line of sight direction. Pulsars, however, 
do not place a specific number of their pulses included in their signals, so each 
pulse 14 can appear to be the same as any other. Pulse cycles from multiple 

52 



pulsars, however, can be combined to help identify which integer cycles best 
match a chosen criteria. The set of integer cycles that most appropriately 
determines the absolute position of the vehicle is selected. 

This process is very similar to the GPS integer cycle ambiguity resolution 
process. Offering an advantage over GPS, pulsars 16 can provide many different 
cycle, or period, lengths. Some cycles are very small (a few milliseconds) and 
others are very large (many thousands of seconds), which assists the pulsar cycle 
resolution method. 
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XII. Velocity Determination Method 

Various applications may require knowledge of a vehicle's velocity, or 
speed and direction. For instance, when a spacecraft requires an orbital 
maneuver, the velocity of the vehicle is used to determine the appropriate point 
in its orbit to accomplish the rocket firing or thrust maneuver. 

To determine the velocity of spacecraft 12, a sequence of position 
estimates can be accumulated using any of the methods described above using 
pulsars 16. The velocity can be computed using the difference of successive 
position estimates divided by the time interval between estimates. Thus, velocity 
represents the differential of position over time. 

Alternatively, the time differentiation of Eqs. (4) and (21) can be 
combined with the time differentiation of Eqs. (25) and (32), to produce an 
analytical expression of velocity based upon pulsar signal measurements. 
Because pulsars 16 transmit pulse signals 14 that are periodic in nature, as a 
spacecraft 12 moves toward or away from the source, Doppler effects will be 
present in the pulse signal 14 due to the spacecraft's motion. For use with 
pulsars, second-order and higher Doppler effects may be present which must be 
included as additional effects in the time differentiation of Eq. (21). 
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XIII. Attitude Determination Method 

Additionally, the attitude of spacecraft 12 may also be determined using 
the observations taken from the pulsar. Assuming a static or fixed detector 22 on 
spacecraft 12, the attitude of the vehicle can be determined by detecting the pulse 
signal 14 and comparing the resulting signal against the database of known 
pulsar profiles. A detector pointed randomly in the sky will either detect a 
recognizable source or the X-ray background. For detected signals above the 
background level, comparisons can be made which help determine which, if any, 
pulsar is in view of the detector. For a static detector that has to wait until a 
pulsar enters its field of view during the orbit of the spacecraft, this process may 
take some time to detect a source. Once a pulsar is detected and identified, 
however, the pulsar location information provides a means of determining the 
sensor, and hence the vehicle, attitude. The line-of-sight to the pulsar will be 
known in inertial frame coordinates, and once detected by the sensor, the sensor 
to inertial frame transformation provides vehicle attitude. Slowly rotating 
spacecraft could also use this type of attitude determination process. It is also 
possible to use X-ray sources that do not pulse successfully within the attitude 
determination process, as long as these sources are uniquely identifiable, and 
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their position coordinates are known. 

Alternatively, a gimbaled sensor system, which allows the sensor to be 
repositioned to point to any point in the sky, can be used to scan various pulsar 
locations in the sky to hasten the process of detecting a signal pulsar. Hanson, in 
his Ph.D. dissertation, "Principles of X-Ray Navigation", 1996, discusses a 
scheme of using gimbaled sensors to determine attitude using X-ray sources. 

XIV. Differential and Relative Navigation Methods 

The use of X-ray pulsars, and other variable X-ray sources, is not limited 
to single spacecraft navigation. For example, an orbiting base station may be 
used to detect pulsar signals and to transmit pulse arrival times and signal errors 
to other spacecraft. This base station can also be used to monitor and update 
pulsar ephemeris information. Ideal locations for these base stations include 
geosynchronous orbits, and Sun-Earth and Earth-Moon Lagrange points. 
Receiving spacecraft are able to navigate with improved information relative to 
the base station. Relative positioning from a main vehicle may also be 
implemented in a satellite formation-flying concept. 

Fig. 1 1 is a schematic view of the alternative embodiment of the method 
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utilizing sources of pulsed celestial radiation for performing differential and 
relative navigation. Similar to Fig. 1 1, a satellite or other spacecraft 100 is in 
orbit of the Earth 20. The satellite or spacecraft 100 has a X-ray detector or 
sensor mounted thereon as in the embodiment of Fig. 1, and is sensing the pulsed 
celestial radiation 14 produced by pulsar 16. In contradistinction to the preferred 
embodiment of Fig. 1, the embodiment of Fig. 1 1 includes a second satellite or 
spacecraft 110 which also receives and senses the celestial radiation 14 produced 
by pulsar 16. In this system, the position of spacecraft or satellite 100 is 
determined with respect to the reference station 110. 

In this system and method of the embodiment of Fig. 1 1 , the incoming 
pulsar phase plane is determined at the reference station 110. The reference 
station then broadcasts the pulsar phase information and, simultaneously, the 
incoming pulsar phase plane is detected by the remote vehicle 100. The phase 
solutions generated by both the reference station 110 and the remote vehicle 100 
are compared in order to determine the relative position of the remote vehicle 
100 with respect to the reference station 110. This differential and relative 
navigation method allows for the effective transmission of navigational 
information to spacecraft in a high-Earth orbit in a manner similar to that of the 
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GPS. Utilizing a system such as this allows for easy expansion to interplanetary 
orbital navigation. 

XV. GPS Complement Navigation Method 

The pulsar-based system can also be used to complement, or provide an 
alternative, to vehicles that are able to use the GPS (and other human-made 
global navigation systems). GPS provides a high quality time signal along with 
transmitting source position information to allow a receiver to precisely 
determine its navigation solution about the Earth. Instances may exist where 
GPS satellites are obscured or unavailable, thus having the pulsar-based 
navigation system as a back-up to critical systems is a benefit. The pulsar-based 
system could further be used as a time-only reference system, which may aid 
high-data rate communication. 

XVI. Sensor Instrument 

X-ray sensors, such as pulse sensor 22 shown in Fig. 6, mounted on 
spacecraft 12, are well-known and have successfully flown on many orbital 
missions, as listed in Table 1 . Although the sensor 22 shown in Fig. 6 is 
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illustrated as being mounted on a satellite, this is for exemplary purposes only 
and the sensor 22 may also be mounted on a spacecraft, planetary rover, or any 
other sort of well-known mobile platform. Various detectors, such as the one 
used by the USA Experiment, are gas-filled collimators used to sense the arrival 
of X-ray photons. Newer sensor technology, such as those based on silicon or 
germanium, can be used as detectors situated at the base of a collimated 
container. A coded-aperture mask or focusing X-ray optics may be used to help 
image the X-ray sources within the field of view. Further, multiple detectors, 
either within the same sensor unit or positioned at strategic locations throughout 
the spacecraft surface, could be used to detect the multiple sources over the same 
time epoch. 

Fig. 8a shows a side view of an X-ray detector 22 which may be mounted 
on a spacecraft 12. The unit is shown as consisting of five detection sub-units, 
however, any number of individual sub-units or sub-sensors may be utilized. 
Each detection unit may include a coded-aperture mask 26. Fig. 8b shows a top 
view of the detector 22 and illustrates the individual coded-masks 26. 

Additionally, each detection unit may include a containment structure, a 
thin collimator, and a silicon-strip detector 28, shown in the bottom view of Fig. 
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8c, positioned directly beneath the collimator. The five detection units are 
positioned with one in a Zenith position and four positioned around this unit at 
45° angles to the Zenith position. The containment boxes shown are 
approximately 10 cm high, 10 cm wide, and 30 cm long, which allows for a 100 
cm 2 detection area. The containment chambers are positioned for reduced 
overall size and the photon paths within each chamber intersect, but no X-ray 
photons or intersecting paths will collide or interact with one another. 

Fig. 6 shows a schematic view of the spacecraft or vehicle 12. A pulse 
sensor 22 is mounted on the spacecraft 12 for detecting X-ray pulses 14 coming 
from a celestial object. The X-ray detector 22 may include a collimator for 
focusing the electromagnetic radiation onto an electromagnetic sensor, such as a 
silicon or germanium sensor. The pulse sensor 22 is in electrical communication 
with a processing means 30 and a memory 34. Processing means 30 is further in 
electrical communication with clock or timer 24. The clock or timer 24 is used 
for measuring the pulse arrival or pulse time period or duration detected by the 
pulse sensor 22. The processing means 30 calculates the navigational 
coordinates of vehicle or satellite 12 utilizing both the timing data and the known 
position of the X-ray sources stored in an X-ray source almanac in memory 34. 
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Power is supplied to spacecraft 12 and the processing means and various 
sensors and control systems by an on-board power supply 32, which may be in 
the form of a storage battery, solar panels, or other power conversion electronic 
system. A vehicle control system 40 is mounted on the spacecraft or vehicle 12 
and is responsive to signals sent by the processing means or through external 
telemetry data, and allows for positioning and directional/motion control of the 
vehicle 12. 

External sensors 38 may also be mounted on the spacecraft or vehicle 12, 
depending upon the specific requirements of the navigational system, and these 
external sensors may include gyroscopes, accelerometers, a GPS 
receiver/antenna, a star camera/tracker, a magnetometer, Sun/Earth sensors or the 
like. The navigation data computed by the external sensors 38 may be used by 
the processing means 30 to further enhance navigation performance. Data from 
external sensors 38 may be blended with data from the pulse sensor 22 and clock 
24 to produce navigation solutions with improved integrity and accuracy. This 
data blending can be accomplished using a Kalman filter, or similar data 
blending algorithm, and can be used in conjunction with the pulse navigation 
flow diagram of Fig. 9. 
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Communications with ground monitoring stations, a base station or with 
other systems may be carried out through communications link 36, which may be 
a radio receiver/transmitter. Telemetry and navigational data may be transmitted 
to and from the spacecraft 12, and also the navigational data calculated by 
processor means 30 may be transmitted to a base station or to other navigational 
systems. 

Fig. 9 is a flow diagram illustrating the computer program aided 
calculation process performed by the navigational system 10. The X-ray photons 
of pulse 14 are received by the detector 22, mounted on the spacecraft 12 of Fig. 
6. The pulse detection signals generated by the pulse signal 22 are sent to 
processing means 30 at which point the calculation process starts at step 42, 
shown in Fig. 9. At step 44, the pulsar's pulse time of arrival (TO A) is 
determined. The pulse time of arrival at the spacecraft 12 is then transformed, or 
calculated, to the time of arrival at the chosen inertial reference frame, chosen to 
betheSSBC in Fig. 1. 

The time of arrival at the SSBC is then compared at step 48 with a pre-set 
data set of a pulsar time of arrival model 52. At step 48, the time of arrival offset 
time is calculated and the navigation solution error is calculated at step 50. If the 
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navigation solution error calculated at step 50 is not within an acceptable range, 
the position, velocity, attitude and time offset data is transferred to step 54 where 
the navigational solution estimate is transferred back to step 46 for a 
recalculation of the time of arrival at the SSBC. This process loop can continue 
with multiple measurements from the same pulsar 1 6, as well as with 
measurements from other pulsars. If the solution error at step 50 is within an 
acceptable range, the system passes to step 56 where either new navigational data 
can be processed by the system, or the navigational process is finished at step 58, 
depending upon the navigational requirements of the user. 

Although this invention has been described in connection with specific 
forms and embodiments thereof, it will be appreciated that various modifications 
other than those discussed above may be resorted to without departing from the 
spirit or scope of the invention. For example, functionally equivalent elements 
may be substituted for those specifically shown and described, and various 
proportional dimensions and quantities of the elements shown and described may 
be varied, and in the method steps described, particular steps may be reversed or 
interposed, all without departing from the spirit or scope of the invention as 
defined in the appended Claims. 
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